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In an external magnetic field B, the spins of the electron and hole will precess in effective fields 
b e + B and b h + B, where b c and b h are random hyperfine fields acting on the electron and hole, 
respectively. For sparse "soft" pairs the magnitudes of these effective fields coincide. The dynamics of 
precession for these pairs acquires a slow component, which leads to a slowing down of recombination. 
We study the effect of soft pairs on organic magnetoresistance, where slow recombination translates 
into blocking of the passage of current. It appears that when b F and b h have identical gaussian 
distributions the contribution of soft pairs to the current does not depend on B. Amazingly, small 
inequivalence in the rms values of &„ and b h gives rise to a magnetic field response, and it becomes 
progressively stronger as the inequivalence increases. We find the expression for this response by 
performing the averaging over fo e , b h analytically. Another source of magnetic field response in the 
regime when current is dominated by soft pairs is inequivalence of the p-factors of the pair partners. 
Our analytical calculation indicates that for this mechanism the response has an opposite sign. 

PACS numbers: 73.50.-h, 75.47.-m 



I. INTRODUCTION 

Due to complex structure of organic semiconductors 
and their spatial inhomogeneity it is nearly impossible 
to identify a unique scenario of current passage through 
them. In view of this, it is remarkable that sizable 
change of current through a device based on organic semi- 
conductor takes place in weak external magnetic fields. 
This effect, called organic magnetoresistance (OMAR), 
seems to be robust, i.e. weakly sensitive to the device 
parameters. Although the first reports on the obser- 
vation of organic magnetoresistance (OMAR) appeared 
decades agoi^, systematic experimental study of this ef- 
fect started relatively recently^— (see also the review 
Ref. 

On the theory side, it is now commonly accepted that 
the origin of OMAR lies in random hyperfine fields cre- 
ated by nuclei surrounding the carriers (polarons) . More 
specifically, the basic unit responsible for OMAR is a pair 
of sites hosting carriers (polarons); the spin state of the 
pair is described by the Hamiltonian 



h = rii • s c + n 2 ■ s h . 



(1) 



Here S a are S h are the spin operators of the pair-partners 
(we will assume that they are electron and hole, respec- 
tively); fli = B + b e and f?2 = B + b h are the full fields 
acting on the spins. They represent the sums of exter- 
nal, B, and respective hyperfine fields, b„ and b h . As was 
first pointed out by Schulten and Wolynes^, due to the 
large number of nuclei surrounding each pair-partner and 
their slow dynamics, b e and b h can be viewed as classical 
random fields with gaussian distributions. 

In order to give rise to OMAR the Hamiltonian 
Eq. ((T|) is not sufficient. It should be complemented 
by some mechanism through which the pair-partners 
"know" about each other, so their motion is correlated 
without direct interaction. The simplest example of such 




FIG. 1: Preferential coordinate system used for analysis of 
the dynamics of the spin pair. Both fields Jli, SI2 reside in 
the rcz-plane. The direction of the quantization axis, z, is 
fixed by the condition fli. x = — fl2,x- 



a mechanism is spin-dependent recombination, i.e. the 
requirement that electron and hole can recombine only if 
their spins are in the singlet, S, state. Then the essence of 
OMAR can be crudely understood as a redistribution of 
portions of singlets and triplets upon increasing B. This 
redistribution affects the net recombination rate. Clearly, 
the characteristic B for this redistribution is ~ b sl b h . 

Naturally, the specific relation between the current and 
recombination rate involves also the rate at which the 
pairs are created. It is important, though, that the latter 
process is not spin-selective. 

Existing theories of OMAR can be divided into two 
groups which we will call "steady-state" and "dynami- 
cal" . The theories of the first group^i appeared earlier. 
In a nutshell (see Ref. [HI for details), in these theo- 
ries the right-hand-side of the equation of motion for the 
density matrix ip = [H , p] with Hamiltonian Eq. ([T]) 
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is complemented with "source" and spin- selective "sink" 
terms. After that, p is set to zero. In Refs. [TH cur- 
rent is expressed via the steady-state p and subsequently 
averaged numerically over realizations of hyperfmc fields. 

The "steady-state" approach applies when the pair 
does not perform many beatings between S and T dur- 
ing its lifetime, since the beating dynamics is excluded 
by setting p = 0. 

This beating dynamics has been incorporated into the 
OMAR theory Ref. which appeared last year. This 
theory relies on decades old findings in the field of dy- 
namic spin-chemistry^. Below we briefly summarize 
these findings. 

If an isolated pair is initially in S, it was shown in Ref. 
that the averaged probability to find it in T after time 
t is given either by the function 



Ps T {t) 



1 
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for strong fields B 3> b c , b h , or by 



Ps T (t) 
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(3) 



for B -C b ei b h . Here b e , b h are the rms hyperfmc fields for 
electron and hole. Naturally, the probability approaches 
3/4 at small B and 1/2 at large B. 

In the theory of Ref. [22| the B-dependent dynamics de- 
scribed by Eqs. ^ translates into the independent 
resistance (OMAR) on the basis of the following reason- 
ing. The dynamics p ST (t) leads to prolongation of the 
recombination time (hopping time, r h , in the language of 
Ref. This prolongation is quantified by 
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(4) 



The meaning of Eq. (gj is that a pair should stay in S 
in order for a hop to take place. Prolongation of hopping 
time leads to a independent increase of the resistance. 
The authors of Ref. [22| evaluated p S r{t) for arbitrary B, 
while in calculation of OMAR they assumed that bare 
hopping times, T h , have an exponentially broad distribu- 
tion. 

Both theories Refs. [l||, [22| take as a starting point a 
pair with the Hamiltonian Eq. ([T]) describing its spin 
states and preferential recombination (hopping) from S. 
The dynamics of this seemingly simple entity, which is 
crucial for OMAR, possesses some nontrivial regimes. 
Uncovering these regimes is a central goal for the present 
paper. The other goal is to demonstrate that nontrivial 
dynamics can manifest itself in OMAR. 

To underline that the spin dynamics of two carriers in 
non-collinear magnetic fields which can recombine only 
from S can be highly nontrivial, we note that separa- 
tion of this dynamics into S-T "beating" stage followed 



by instantaneous hopping after time T h , as in theory |22j, 
is not always possible. It is quite nontrivial that spin- 
selective recombination of carriers can exert a feedback 
on the spin dynamics. As an illustration of this deli- 
cate issue we invoke the example of cooperative photon 
emission discovered by R. H. Dicke^. In the Dicke ef- 
fect one superradiant state of a group of emitters having 
a very short lifetime automatically implies that all the 
remaining states are subradiant and have anomalously 
long radiation times. Below we demonstrate that a sim- 
ilar situation is realized in dynamics of two spins when 
recombination from S is very fast. We will see that the 
remaining 3 modes of the collective spin motion become 
very "slow" . 

Our analysis reveals the exceptional role of the "soft" 
pairs, which are sparse configurations of b el b h for which 
full fields S~2i , 1^2 have the same magnitude. 

The paper is organized as follows. In Sect. II we cast 
the eigenmodes of the Hamiltonian Eq. ([1]) in a conve- 
nient notation. In Sect. Ill we include recombination and 
study its effect on the eigenmodes. The consequences of 
nontrivial dynamics for OMAR are considered in Sects. 
IV and V, where we perform averaging over realizations 
of hyperfmc fields. We establish that inequivalence of rms 
hyperfine fields for electrons and holes has a dramatic ef- 
fect on OMAR, when it is governed by soft pairs. Sect. 
VI concludes the paper. 



II. DYNAMICS OF A PAIR IN THE PRESENCE 
OF RECOMBINATION 

A. Isolated pair 

We start with reviewing the dynamics of a pair of spins 
in the absence of recombination. Obviously, this dynam- 
ics does not depend on the choice of the quantization 
axis. However, since we plan to include recombination, 
the choice of the quantization axis, z, illustrated in Fig. 
Q] appears to be preferential. The axis is chosen to lie 
in the plane containing the vectors fli, Q,2- Moreover, 
the orientation of the z-axis is fixed by the condition 
tti x = — Q.2X- Then the angles, 6\, 82, between tt 2 
and the z-axis are given by 



tan 6*i 



tan( 



(5) 



With this choice, the Schrodinger equation for the am- 
plitudes of S, Tq, T+, and T_ reduces to the system 



.0T+ 1 
dt + y/2 

i-rr- = A^Tq - -^A X T+ + — 7= Aj;T_ 
dt ° V2 V2 



.dT 



(6) 
(7) 
(8) 
(9) 



3 



where S 2 , A z , and A x arc denned as 
A, 
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The advantage of our choice of the quantization axis 
shows in the fact that the state To is coupled exclusively 
to S. Since recombination is allowed only from S, this 
will simplify the subsequent analysis of the recombina- 
tion dynamics. 

The eigenvalues, Aj, of the system Eqs. ©-© satisfy 
the quartic equation 



A?(A?-E2)-A?(A 



0. 



(13) 



We will enumerate these eigenvalues according to the 
convention Ai = — A2 and A3 = — A4. To find the absolute 
values Af, A3 one does not have to solve Eq. (|13p. since it 
is obvious that for non-interacting spins the eigenvalues 
are the sums and the differences of individual Zeeman 
energies 



A? = 



ini| + in a 



\ 2 - 
a 3 - 



inj- in a | 



(14) 



Naturally, Ai, A3 do not depend on the choice of the 
quantization axis. At the same time, the coefficients in 
Eq. ([13]) do depend on this choice. To trace how the 
dependence on the quantization axis disappears in the 
roots of Eq. (fl~3|) . one should use the following identities 



A 2 = Iftll 



2 

l^il 2 - 



(15) 
(16) 



In terms of the angles 0\ and #2, Fig. [TJ the corre- 
sponding eigenvectors can be expressed as 



condition we also have d\ — 82- Then the first two eigen- 
vectors Eq. (|17p have zeros in the rows corresponding to 
To. Concerning the other two eigenvectors, due to their 
degeneracy, their sum and difference are also eigenvec- 
tors. The difference has a zero in the row corresponding 
to T), while the sum consists of the T) component, exclu- 
sively. Then we conclude that for realizations of hyper- 
fine field for which = \ fl2\ the state Tq is completely 
decoupled from the other three states. This fact has im- 
portant implications for recombination dynamics, as we 
will see below. 

Including recombination requires the analysis of the 
full equation for the density matrix 



ip=[H,p]- — {p,\S) (S\}, 



(18) 



where r is the recombination time. The form of the sec- 
ond term ensures that recombination takes place only 
from S. The matrix corresponding to Eq. (fT8]l is 16 x 16. 



The 16 eigenvalues can be cast in the form A; — A* , where 
Ai and A, satisfy the equation 



Aj A, + - (\j -Hi)- Xf(A z + Ai) + Ai^ z = 0. (19) 



The latter equation expresses the condition that Ai are 
the eigenvalues of non-hermitian operator H ~ ^\S) (S\. 
In the limit t — > 00 this equation reduces to Eq. (| 14[) . 
The dynamics of recombination is governed by the imagi- 
nary parts of the roots of Eq. (fl"9|). i.e. decay is described 
by the exponents exp [— (Im Ai + Im Xj ) t] . Less trivial 
is that finite r can strongly affect the real parts of A;. 
Physically, the dependence of ReAi on r describes the 
back-action of recombination on the dynamics of beating 
between different eigenstates. In the following two sub- 
sections this effect will be analyzed in detail in the two 
limiting cases. 
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where the first two correspond to Ai,2 while the last two 
correspond to A34, respectively. 

The form Eq. (fTT)) allows us to make the following ob- 
servation. When the full magnetic fields acting on spins 
incidentally coincide, we have = | 1 . Then it fol- 
lows from Eq. (fT4"|) that A3 = A4 = 0, so that the two cor- 
responding eigenstates become degenerate. Under this 



B. Slow Recombination 

Consider the limit 1/r < A ; . In this limit recombina- 
tion amounts to the small corrections to the bare values 
of Ai given by Eq. (fl"4"]) . This allows one to set Ai equal 
to their bare values in all terms in Eq. (|19p containing 
1/r, and search for solution in the form A; + <5A.;. Then 
one gets the following expression for the correction SXi 



SX,, = — 



Ai(Af-S|) 



-i Af(Af-E^) 

r Af-2A l (S2 + A 2 + A2) 2r Xj - J%A* ' 

(20) 

In the last identity we have used the fact that A; satisfy 
the equation Eq. (|14l) . The above expression can be 
greatly simplified with the help of the relations Eq. (|15|) . 
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FIG. 2: (Color online). (a) Slow-recombination regime, 
fii,2 <C 1/t. Horizontal lines represent the energy levels Eq. 
(14) of a pair in the absence of recombination. Recombina- 
tion from S causes the broadening of the levels Eqs. (21,22), 
which, for a typical pair, is of the same order for all levels, (b) 
Slow-recombination regime. For soft pairs, |fii| ~ IO2I, re- 
combination results in splitting Eq. (24) of the widths of the 
levels A3,4 rather than their positions, (c) Fast-recombination 
regime, S7i,2 3> 1/t. The eigenstates S, To, T + , and T_ are 
well-defined. Recombination causes strong broadening, 1/t, 
of the level S, and weak broadening ~ fii,2T of the other three 
levels. 



One has 
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^34 
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(21) 
(22) 



The above result suggests that for generic mutual ori- 
entations of fix and O2 all modes of a pair decay with 
characteristic time ~ r. At the same time, for parallel 
orientations of fix, SI2 the modes A12 have anomalously 
long lifetime. This long lifetime has its origin in the fact 
that for fix || fi2, the states T+ and T_, which arc orthog- 
onal to S, are the eigenstates of the Hamiltonian Eq. (QJ . 
Formally this can be seen from the general expression Eq. 
(fTT| for the eigenvectors upon setting 0\ = 62 =0. Sim- 
ilarly, for fix and fi2 being antiparallel, one can check 
from Eq. (|17[) that for 6\ = tt — #2 , the eigenstates A3 , A4 
have no S component, so they are long-lived. Note that 
the existence of long lifetimes for parallel and antiparallel 
configurations of fix, fi2 is at the core of the "blocking 
mechanism" of OMAR proposed in Ref. [l5l 



1. Soft pairs 

As was pointed out in the Introduction, recombination 
also has a pronounced effect on the spin dynamics for 
sparse configurations for which |f2x| ~ |fi 2 |- Indeed, for 
these configurations, the values A3 and A4 are anoma- 
lously small. Then the basic condition, 1/t <C Ai, under 
which Eq. (|2"Tj) was derived, is not satisfied. We dub such 
realizations as soft pairs. For soft pairs the expressions 
for SXi, SX2 remain valid, but the eigenvalues A3, A4 get 
strongly modified due to finite recombination time, r. 

Although for soft pairs the terms oc 1/t in Eq. (fH))) 
cannot be treated as a perturbation, a different simplifi- 
cation becomes possible in this case. We can neglect A? 



compared to E 2 , in the first term and A 2 compared to A 2 
in the second term. The first simplification is justified, 
since the typical value of S z is ~ |fix| ~ |fi2| and is much 
bigger than both 1/t and (|fii| — |fi2|)- Concerning the 
second simplification, the smallness of (|fix| — |fi2|) au- 
tomatically implies that A z given by Eq. (fT0|) is small. 
With the above simplifications the eigenvalues A 3i 4 sat- 
isfy the quadratic equation 



(E2 + A£)A? + -E2A i 



2ip2 



A^S 
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(23) 



Already from the form of Eq. (|2"3"|) one can make a sur- 
prising observation that, even with finite 1/t, one of the 
roots is identically zero when A 2 = 0, i.e. when |fix| and 
|fi2| arc exactly equal to each other. This suggests that 
a pair in the state corresponding to this root will never 
recombine. For a small but finite difference (|fii| — |fi 2 |) 
the recombination will eventually take place but only af- 
ter time much longer than r. Indeed, for the generic case, 
I fix I ~ |fi 2 |), we have from Eq. $Z5§ 



A 3 ,4 - — 



A ± v/A 2 - 4AA 2 t 2 



where the dimensionless parameter A is defined as 



A 
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(24) 



(25) 



Even when |fix| and |fi2| are close, a typical value of pa- 
rameter A is ~ 1. Then Eq. (|24]l suggests that anoma- 
lously long- living mode exists in the domain A 2 < 1/t 
where its lifetime is ~ 1/A 2 t. Note that the lifetime be- 
comes longer with a decrease of the recombination time. 

As the difference |fix| — |fi2| increases, the product 
A z t becomes big and the expression under the square 
root in Eq. (|24|) becomes negative. Then the lifetimes of 
of both states corresponding to A3 and A4 become equal 
to t/A. Note that, at the same time, the splitting of the 
real parts of A3 and A4 becomes ~ A 2 t, which is much 
bigger than |fix| — |fi2|- 

The above effect can be interpreted as a repulsion 
of the eigenvalues caused by recombination^. A more 
prominent analogy can be found in optics^. The signs 
+ and — in Eq. can be related to the superradiant 
and subradiant modes of two identical emitters. The role 
of t in this case is played by their radiative lifetime. 

Both effects illustrate the back-action of recombina- 
tion on the dynamics of the pair when the spin levels of 
pair-partners are nearly degenerate. To track an analogy 
to this effect one can refer to Refs. and [2f| where 
Eq. (|24p appeared in connection to resonant tunneling 
through a pair of nearly degenerate levels, while the role 
of 1 /t was played by the level width with respect to es- 
cape into the leads. 

For our choice of the quantization axis the long-living 
state corresponds to To. For completeness we rewrite the 
parameter A z , which enters Eq. (1241) . in the coordinate- 
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independent form 
A 
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To establish coordinate-independent form of parameter 
A we need the combinations £ 2 and £ 2 + A 2 , which are 
given by 



^ 2 _ |«i + n 2 | 
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so that A can be cast into the form 



A 
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(27) 
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The consequences of "trapping" described by Eq. (pM)) 
for OMAR will be considered in Sections IV and V. In 
the subsequent subsection we will see that the similar 
physics, namely, the emergence of slow modes due to fast 
recombination persists also in the domain |J7i j2 |t <C 1. 



C. Fast Recombination 

In the opposite limit, r -c |Jli. 2 | _1 , the bracket (Xi + ^) 
in Eq. (fT9")l is big. This suggests that three zero-order 
eigenvalues are 



Xi =o,±s z . 



(30) 



In the same order, the fourth eigenvalue is — —. Concern- 
ing the eigenvectors, in the zeroth order they are simply 
S, T + ,T_, and To. This follows from the equation 

iS + -S = A Z T - -^T+ + -Lt_. (31) 



Taking t to zero means that in the zeroth order 5 = 
Then three other equations in the system Eq. © get 
decoupled. 

In the first order, the eigenvalues Eq. (|30p acquire 
imaginary parts 
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With the help of Eqs. (|2"T1) and (|2U)) these imaginary 
parts can be simplified to 
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We see that for a generic situation |f2i| ~ |n 2 | the life- 
time of the modes To, T + , and T_ are ~ 1/|J~2i. 2 | 2 t, i.e. 



in the regime of fast recombination it is much longer than 
t. This is a consequence of effective decoupling of To, T + , 
and T- from S in this regime. We also observe from Eq. 
(|3"3")l that there is additional prolongation of lifetime for 
the mode To if the pair is soft. Eq. (|33|) also suggests 
that lifetimes of the states T + , T_ arc anomalously long 
when Jli and f2 2 are collincar. This expresses the obvi- 
ous fact that, for collinear effective fields acting on the 
pair-partners, T + and T_ are the eigenstates no matter 
whether recombination is present or not. 

Once the eigenvalues and eigenvectors of a pair in the 
presence of recombination are established, the next ques- 
tion crucial for transport through the pair is: Suppose 
that initial state is a random superposition of S, T , T + , 
and T_ , what is the average (over the coefficients of su- 
perposition) waiting time for this state to recombine? 
Naturally, the answer to this question does not depend 
on the actual choice of the orthonormal basis. We ad- 
dress this question in the next section. 



III. RECOMBINATION TIME FROM A 
RANDOM INITIAL STATE 

A. Soft pair in a slow recombination regime 

To illustrate the peculiarity of the question posed 
above, we start from an instructive particular case of 
soft pair in a slow recombination regime. We defined a 
soft pair as a pair for which the condition (|J7i| — |r2 2 |) <C 
|fii 2 1 is met. However, in the slow recombination regime, 
the combination (\fti\ — \Sli\)T can be either big or small. 
In both cases there is a strong separation between the ab- 
solute values of Ai j2 and A3 ; 4. It can be seen from Eq. 

that in the limit - |n 2 |)r > 1, the recombina- 
tion times for states which correspond to A3 and A4 are 
given by 



— and 
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get while in the opposite limit, (|Oi| — |J^ 2 |)t <C 1, we get 



2r 
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and t { . 
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(36) 



We see that the recombination time of A3 is ~ r for 
both limits, while the recombination time of A4 crosses 
over from ~ t to ~ 1/A 2 r as — |J~2 2 |)t decreases. 

Taking into account that for generic case ~ |fi 2 | the 
recombination times corresponding to Ai j2 are ~ r, we 
conclude that for purely random initial conditions the 
average recombination time is either ~ r or it is j of 
1/A 2 r. 

The major complication for getting exact average re- 
combination time for a soft pair is that the exact eigen- 
states represent mixtures with weights governed by the 
recombination time. This follows from Eq. (|24[) . In addi- 
tion, the eigenstates corresponding to A3 , and A4 are not 
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orthogonal to each other. However, for a soft pair these 
complications can be overcome. The reason is that, there 
are two small parameters in the problem, l/r|fii. 2 |, and 
( | ill | — | ^2 1 ) / 1 Oi 2 1 • The first parameter guarantees slow 
recombination, while the second ensures that the pair is 
soft. The presence of these parameters allows us to eval- 
uate (t R ) in the closed form using the general formula 



i 



(37) 



where gij = (vi\vj) is a matrix of inner products of eigen- 
vectors corresponding to complex eigenvalues A^ and A,,- . 
The above formula becomes absolutely transparent when 
the eigenvectors are orthonormal. Then the matrix 
reduces to the Kronecker symbol, <5y, and (t R ) simplifies 
to 



1 



Im A, 



(38) 



which expresses the fact that for random initial state the 
average recombination time is the evenly-weighted sum 
of recombination times from eigenstates. 

In the case of a soft pair and slow recombination one 
should use Eq. (|3"T|) to evaluate (t R ). What enables 
this evaluation is that, by virtue of small parameters, 
the eigenvectors corresponding to Ai and A2 are mutu- 
ally orthogonal (with accuracy 1/t| Oi 3 2 1 ) , and they are 
both orthogonal to eigenvectors corresponding to A3 and 
A4. Therefore, in evaluating Eq. (|37p . one has to deal 
only with mutual non-orthogonality of two eigenvectors 
V3 and V4. The straightforward calculation yields 
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A 4A?r ImAi ImA 2 



(39) 



where ImAi = I1T1A2 are given by Eq. (|21[) . It is 
easy to see that in the limiting cases of large and small 
(|rii| - |« 2 |)t Eq. ([33) reproduces Eqs. {35J and ([55]). 
respectively. 

While in the last two terms in Eq. (|39| depend weakly 
on the degree of "softness" of the pair, A z cx — |fi 2 |), 
the second term exhibits unlimited growth with decreas- 
ing A z . We emphasize the peculiarity of this situation. 
In conventional quantum mechanics, when the level sep- 
aration becomes smaller than their width, it should be 
simply replaced by the width. What makes Eq. ((39)) 
special is that the smaller is A z the more the state Tq 
becomes isolated. There is direct analogy of this situ- 
ation with the Dicke effect^, as was mentioned in the 
Introduction. By virtue of this analogy, the state To as- 
sumes the role of the "subradiant" mode which accom- 
panies the formation of the superradiant mode. In the 
Dicke effect the formation of superradiant and subradi- 
ant states occurs because the bare states are coupled via 
continuum. In our situation it is recombination that is re- 
sponsible for "isolation" of T . If the pair is not soft, the 
calculation of the time (t R ) in the slow- hopping regime 
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FIG. 3: (Color online). The simplest model of transport 
through a bipolar device in which the currents flow along inde- 
pendent chains. Electrons arrive at the recombination region 
from the left, while the holes arrive from the right. Blobs en- 
close the sites from which electron and hole recombine. One 
of the blobs is enlarged to illustrate the spin precession of the 
pair partners in their respective fields Oi, fl2- For soft pairs 
the magnitudes of fli and {I2 are close to each other. 



can be performed by simply using Eq. (|3"8|) and Ai given 
by Eqs. pTj) . ([24]). This is because the smallness of 
1/t makes the eigenstates almost orthogonal. However, 
the Dicke physics becomes even more pronounced in the 
fast-recombination regime, as demonstrated in the next 
subsection. 



B. Recombination time in the fast recombination 
regime 

It might seem that under the condition of fast recom- 
bination |S7i, 2 |t *C f the recombination time from the 
random initial state should be ~ r, since spins practically 
do not precess during the time r. The fact that recom- 
bination takes place only from S, while initial state is a 
random mixture, already suggests that (t R ) is longer than 
t. This is because if the initial configuration is different 
from S it must first cross over into S by spin precession 
before it recombines. The characteristic time for the spin 
precession is ~ | $"2 1 .2 1 1 ^ t. It turns out that the cross- 
ing time is actually much longer than |Oi 2I • Formally, 
this fact follows from Eqs. (|3"3"|) . (|34[) for 5 A,, which are 
of the order of |f2i )2 | 2 r rather than | $~2 1,2!- We can now 
interpret this result by identifying S with superradiant 
state, while To, T + , and T_ assume the roles of subra- 
diant states. The short lifetime of S isolates it from the 
rest of the system. Quantitatively, the portion of S in 
the other eigenvectors is ~ |f2i, 2 |T. 

What is important for calculation of (t R ) is the 
fact that eigenvectors are orthogonal (with accuracy ~ 
1 / 1 .2 1 '7") in the fast-recombination regime. This allows 
one to replace the overlap integrals g tj in Eq. (|37p by 5 iS 
and use the Eq. (|38[) which immediately yields for (t R ) 
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FIG. 4: (Color online). I, II, III, and IV are possible variants 
of the current cycle. For each variant the pair is initially cre- 
ated in one of four states. This is followed by time evolution, 
illustrated by blue double arrows, which mixes the states. 
Subsequently, the pair either recombines from S (brown ar- 
row) or dissociates. The processes of creation and dissocia- 
tion are indicated by white double arrows. The current is the 
inverse duration, t, of the cycle averaged over initial states, 
which we assume to have equal probabilities. The time, t, is 
given by Eqs. (44), (45), or (46) depending on the recombi- 
nation regime. 



the result 



(**> 



Im Ac 



ImA 



To 



Im An 



ImA 



T- 



1 



1 

r I A? 



4 

A 5 " 



(40) 
(41) 



Substituting the coordinate-independent expressions for 
A x and A z , we arrive at the final expression for recom- 
bination time, which is applicable within the entire fast- 
recombination regime 



current is due to electron-hole recombination. As shown 
in Fig. [3J electrons arrive at the pair of sites (enlarged 
regions in Fig. |3j from the left, while holes arrive from 
the right. Once an electron-hole pair is formed, the spins 
of the pair-partners undergo precession in the fields SI i 
and f22, respectively, waiting to either recombine or to 
bypass each other and proceed along their respective cur- 
rent paths. For simplicity we choose the current paths in 
the form of ID chains. This choice makes the adopted 
model of transport very close to the "two-site" model 
proposed in Ref. [TEl The on-site dynamics of a pair with 
recombination was studied in detail in previous sections. 
To utilize the results of Sect. Ill for the calculation of 
current, /, one has to incorporate the stages of formation 
and dissociation of pairs into the description of transport. 

In Fig. [4] the formation and dissociation are illustrated 
with white double-sided arrows. The formation time for 
all four variants of initial states is assumed to be the 
same, r D . For simplicity we choose the average time for 
bypassing to be also t d . Note that this choice does not 
limit the generality of the description, provided that r D is 
longer than the recombination time. The middle and the 
bottom portions in Fig. [4] illustrate the spin precession 
(blue arrows) and recombination (brown arrow) stages, 
which we studied earlier. Implicit in Fig. 21 is that the 
pair disappears either due to dissociation or by recom- 
bination before the next charge carrier arrives. Another 
way to express this fact is to state that the passage of 
current proceeds in cycles. 

Naturally, subsequent cycles are statistically indepen- 
dent. This allows one to express the current along a ID 
path through the average duration of the cycle, t. In- 
deed, N ^> 1 cycles take the time T N = t\ + 1% -| h i„ . 

For large N, this net time acquires a gaussian distribu- 
tion centered at T N = Nt. Correspondingly, the current, 
N/T N , saturates at the value 



(t*) 



T + 



lo-i + n. 



\th + ngj 
\Ui x n 2 | 



(42) 

As was already noticed in the previous section, recom- 
bination time diverges for two particular configurations: 
soft pairs with = | O2 1 and collinear and $^2 • 

Certainly this divergence will be cut off in the course of 
calculation of current through a pair to which we now 
turn. 



IV. TRANSPORT MODEL 

We adopt a transport model illustrated in Fig. [3J For 
concreteness we will discuss a bipolar device, so that the 



(43) 



Note, that Eq. (|43|) constitutes an alternative approach 
to solving the system of rate equations for two-site model, 
as in Ref. [lH or to solving numerically the steady-state 
density-matrix equations, as in Ref. \m Note also, that 
Eq. (|43]) is applicable to such singular realizations as soft 
pairs, while previous approaches are not. For detailed 
discussion of this delicate point see Ref. [2ft 

The remaining task is to express t via the average re- 
combination time, (t R ) and r D . For a typical pair in the 
regime of slow recombination (t^ ) is given by Eq. (|3"5]) 
upon substitution of Eq. (|2ip . Using this expression we 
get for average duration of the cycle 



The first term captures the formation of the pair, while 
1/t d in the denominators describes the bypassing. In- 
deed, if recombination times are ~ r, one can neglect 
1/t d in the denominators. On the other hand, as the 



brackets in denominators in Eq. (|44|) turn to zero, which 
corresponds to anomalously slow recombination, the sec- 
ond term becomes t d . Similarly, for slow recombination 
with soft pairs, using Eq. (|39|) we get 



t = T n 



1 
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(45) 



Finally, in the regime of fast recombination one should sion for t 
use Eq. (|42j) for (t R ) . This leads to the following expres- 



t = Tr. 
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1 




1 




4 






l"2| 2 ) 2 \ 

hfl 2 | 2 J 


+ ^7 



2x 



T f |r2 lX r2 2 | 2 \ 




2 \^ |S^i+«2 | y 7 





(46) 



Obviously, the dependence of current on external field 
is encoded in Eqs. (|4"4")) - (|4"6"|) via the frequencies fli = 
B + b e and Q2 = B + fo/j. The observable is the current 
averaged over realizations of the hypcrfine fields b e and 
fo/j. This averaging is performed in the next section. 



V. AVERAGING OVER HYPERFINE FIELDS 

A. Averaging in the slow-recombination regime 

Our basic assumption is that the time, r D , of forma- 
tion and dissociation of a pair is much bigger than the 
recombination time, r. Only under this condition the 
pair will exercise the spin dynamics. Using the relation 
t d ^> t, we can simplify the expression Eq. (|31)l for t of 
a typical pair 



t = r n 



1 - 



|Oi||n 3 | 



(47) 



We can also rewrite the current in the form I = — 

td 

8I t (B), where the field-dependent correction is defined 
as 



Slt(B) = 



1 - 



I fii||n 2 | 



(48) 



As we will see below, the significant change of 6 It with B 
takes place in the domain where B is much bigger than 
the hyperfine field. Therefore, we expand Eq. (|48|) with 
respect to \b e \/B and \b h \/B. The principal ingredient of 
this step is the expansion of denominator 



|ni| 3 |n2| 2 -(ni.n 2 ) s 



B 



Assuming identical Gaussian distributions of 6 C , b h 
1 



■expHWo), 



and choosing the z-direction along B we get 

1 



(49) 



(50) 



(Sh(B)) = 



B 2 t 



{bi x -b 2x f + (b ly -b 2y ) 2 + ^-B 2 



(51) 
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FIG. 5: (Color online). Blue line: Magnetic field response, 
Slt(B), for the "parallel-antiparallel" blocking mechanism, is 
plotted from Eq. (|54[) in the units 1/t d versus dimension- 
less magnetic field B/B c . Green line: fit with conventional 
lineshape of OMAR, x 2 /(0.S + x 2 ). 



The next step is averaging Eq. (1511) over the remaining 
four components of the hypcrfine fields. It is easiest to 
perform this integration by switching to b\ ± bi and in- 
troducing the polar coordinates. The integrations over 
the sum and over the polar angle are elementary. The 
result can be cast in the form 



(SI t (B)) = -jYf- 
t d \B € 



where the characteristic field B c is given by 

1/2 



B a = 



2t d 

T 



b . 



(52) 



(53) 



On the physical level, the fact that the "body" of magne- 
toresistance lies in the domain B ^ b suggests that the 
origin of the effect are trapping configurations for which 
f2i and SI2 are almost parallel or antiparallel. In this 
regard, Eqs. (|52|) and (|54| can be viewed as analytical, 
rather than numerical, as in Ref. [TEl treatment of the 
bipolaron mechanism^. 



B. Averaging in the soft-pair-dominant regime 

Soft pairs are responsible for the second and third 
terms in the brackets of Eq. (|45[) for t . The second term 
becomes big when the sum, fli + Q2, becomes anoma- 
lously small. Still it cannot dominate over the contribu- 
tion from the first term for the following reason. When 
f2i + ^2 is small, the expression in the parenthesis of the 
second term behaves as (fii + 2 ) 2 / 1 Oi | 2 . At the same 
time, for small f2i — H2, the expression in the paren- 
thesis of the first term behaves as (f2i — O2 ) 2 / 1 1 1 2 - In 
strong fields, the second expression is smaller than the 
first, leading to the larger SI, while in weak fields the 
two expressions give the same contribution to SI. 

The third term in Eq. (|45[) captures the contribu- 
tion of the slow modes to the current. Below we will 
study whether the averaging of this term over hypcrfine 
fields can dominate over the "bipolaron" magnetic-field 
response given by Eq. (|55|) . 

Prior to performing averaging, we rewrite the current 
as I = — 5I S (B), like we did above. In the soft- 
pairs-dominatcd regime the expression for SI S (B) takes 
the form 



SIs(B) 



The form of the function J- is the following 



Tn 



|Oi+n 2 | 2 



(56) 



F{x) = 2x 2 



du ■ 



x 2 e x 



Ei(x 2 ), (54) 



where Ei(z) is the exponential integral function. From 
Eq. (|5"5|) we see that relation t d ^> r ensures that B c >• 
bo, so that the expansion Eq. (l49l) of Slt(B) with respect 
to hyperfine fields is justified. 

The magnetoresistancc Eq. (|52[) is plotted in Fig. 
[5j We note that the shape, being a single-parameter 
function, T{x), can be very closely approximated with 
x 2 /(0.8 + x 2 ). This approximation, which is also plotted 
in Fig. [5j represents a standard fitting function for ex- 
perimentally measured magnetoresistancc. It can be seen 
that at x <C 1 there is a small deviation of T(x) from the 
approximation. This is due to singular behavior of J-(x) 
at small arguments. This singularity translates into the 
following behavior of Slt (B) 



Slt(B) 



b <B<B c 



TB 2 
O t 2 U2 

LT D °0 



hi 



2r D b 2 
tB 2 



(55) 



For a typical configuration with | $~2 x | ~ |^2|, the second 
term in denominator can be estimates as Ifiil 2 



'TTr 



that it is large in the slow-recombination regime. This is 
why the soft pairs with 



|«i| - |n 2 |) 



1 



(57) 



give the major contribution to the average SI S (B). The 
latter fact allows one to simplify the averaging proce- 
dure. Namely, one can use the fact that for e -^i 1 the 
combination e 'i+ x z can be replaced by nS(x). Thus, the 
expression to be averaged can be rewritten in the form 



SI S (B) 



|fti| 2 H» 2 | 2 

|Oi +f2 2 l 



(58) 



The form Eq. (|58|) suggests that characteristic magnetic 
field determined from zero of the S- function is B ~ 60, 

and yields the estimate I/t 1 / 2 -^ 2 ^ for SI S (B). To com- 
pare the contribution of soft pairs to that of typical pairs 
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FIG. 6: (Color online). Magnetic field response for the "soft- 
pair" mechanism is plotted from Eq. Q66p versus magnetic 
field in the units of the hyperfine field i>i for different values 
of the asymmetry parameter r\. Inset: fit of the response in 
the limit of strong asymmetry with conventional lineshape of 
OMAR, 72x7(0-23 + a; 2 ). 



The six-fold integral Eq. (|60|) can be reduced to a 
single integral in three steps. As a first step, we introduce 
new variables v = b c — b h and u = b c + b h + 2B, so that 
Eq. (f6"0"|) acquires the form 



(SIs(B)) 



d 3 u J d 3 v \u\ S(u ■ v) 



8ir 2 b 3 b 3 2V ^ 
x exp (-a(u - 2B) 2 + f3(u - 2B) ■ v — a\v\ ) , (61) 
with parameters a and (3 defined as 



1 ( 1 



4 U? + bl ) ' 13 ~ 2 U bl 



(62) 



As a second step, we perform integration over the vector 
v. The reason why this integration can be carried out 
analytically is that, upon choosing the z-direction along 
u, the (5-function fixes v z to be zero. The remaining two 
integrals over v x and v v axe simply gaussian integrals, so 
we get 



this estimate should be compared to Eq. (|55|) taken at 
B ~ bo- Soft pairs dominate if the condition 



— > b T 

T 



(59) 



is met. Since t d is much bigger than r, this condition 
is compatible with the condition, b$r 3> 1 necessary for 
slow recombination. Note in passing, that replacement 
of the denominator in Eq. (|48p by a (5-function, as we did 
for soft pairs, is not permissible. This follows, e.g., from 
Eq. (|53f which suggests that the characteristic field B a 
is much bigger than bo. Replacement of the denominator 
in Eq. (|48|) by a J-function would automatically fix the 
characteristic field at B ~ &o- 

In averaging of Eq. (|58|l over hyperfine configurations, 
we will assume from the outset that the characteristic 
hyperfine fields, b\ and 62, for the electron and hole are 
different, so that 



(SUB)) = 




^b\b\ 
£?| 2 - + i?| 2 
+ b h +2B\ 2 



d 3 6 c / d 3 b h 



exp 



Ir 
"2 



(60) 



Subsequent analysis will indicate that different b\ and 62 
is a necessary condition for 8I S to exhibit independence. 



(SIs(B)) = 



exp 



Sirbfb^a^ tt 

2 , 



d 3 u 



-a(u-2B 




(B-u) 



u 



(63) 



To perform the integration over it, we switch to spherical 
coordinates with polar axis along B. Then the integra- 
tion over azimuthal angle reduces to multiplication by 
27T. The third step is the integration over the polar angle 
in Eq. (J63]) . We have 



(SI S (B)} 



Ab^ay/Trl 



du u e 



2„-au /L 









d9 sin 9 exp ( AauB cos 9 B 2 cos 2 1 

a 



(64) 



Now we note that the integral over 9 can be expressed 
via the error-functions in the following way 



dx e 



V 7T c 

2A 



erf [A+ 



C 



erf A 
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(65) 

We are left with a single integral over u, which can be 
cast in the form 



(SIs(B)) = 



8b 3 b 3 , 



r 3 V^0B 



du u 2 exp ( —a ( 1 --r- 
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(66) 
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The result Eq. (|66|) can be recast in the more con- 
cise form in terms of the Dawson function D(x) = 

" X 2 
t 2 



e~ x Jdte 
o 



The corresponding expression reads 



(6I a (r],B)) V2 



n 



V2 

(51.(0,0)) y/Z=j [2-r,J 2z 
where we have introduced z = B jb\. 



D 



2- 



. (71) 



In the limit of strong asymmetry, when r\ is close to 1, 
one gets a simple analytical expression for (8I S (B)) 



FIG. 7: (Color online) Magnetic field response caused by the 
difference, in the (/-factors of electron and hole is plotted from 
Eq. (|73[) for several values of relative difference, k. Upper in- 
set illustrates that the shape of the response is near-gaussian. 
Lower inset illustrates that at « close to 1 the shape of the 
response develops a maximum. 



(6I 3 ( V = 1,B)) 

(si.(v = o)) 



21 if) I dxvT 

-1 



x exp 



(72) 



C. Analysis of Eq. (|66|) 

At this point we make an observation that for b\ = b 2 , 
which is equivalent to f3 = 0, magnetic field drops out of 
Eq. ([55)) . The easiest way to see it is to set /3 = at the 
earlier stage of calculation, namely in Eq. (|63p 



(SI s (h=b 2 )) 




d 3 ue - a ( u -2B) 



(67) 

which is clearly independent of B after a simple coordi- 
nate shift. If we set b\ = b 2 , then 5I S is given by 



(SI s (h =b 2 )) 



1 



2T7-3 bl ' 



(68) 



in agreement with the qualitative estimate above. 

Magnetic field dependence of (SI S ) emerges already at 
small values of asymmetry parameter defined as 



'/ 



1 



bj 
bl 



(69) 



This is illustrated in Fig. [5J where (SI S (r), B)) — (SI S (rj, 0)) 
in the units of (5I s (r) = 0)), given by Eq. (|rJ8"|) . is plot- 
ted for several values of r]. We see that, as r\ increases, 
the shape of the curves does not change much. For the 
saturation value the analysis of Eq. (jof)]) yields 



(6I s (r),oc)) - (6I 3 (ri,0)) V2 V 2 



(si s (v = o)) 



(2 - t?) 5 / 2 



(70) 



At small B the ratio Eq. ([72]) behaves quadratically, 
while at large B it saturates as \/2 ^1 — g^i^ • Over- 
all, similarly to It(B), magnetoresistance Eq. (|72[) can 
be closely approximated with v / 2^ 2 /(0.23 + x 2 ), as illus- 
trated in Fig. [51 



D. Inequivalence of electron and hole g- factors 

In the previous subsection we demonstrated that ex- 
ternal magnetic field drops out from the general expres- 
sion Eq. (|60[) when the variances b\ and b 2 are equal. 
Here we note that averaging does not eliminate the in- 
dependence even when b\ = b 2 , as long as the g- factors 
of the pair partners are different. Incorporating g\ and 
g 2 into Eq. (|60j) is straightforward and amounts to mul- 
tiplying b e + B by 1 + k, while b h + B is multiplied by 
1 — ft, where k is the relative difference in the g-factors. 
The three steps leading from Eq. (joT))) to Eq. (jfJrJ)) are 
exactly the same as for k = 0. Finite k modifies both 
the prcfactor in the integral Eq. (|66[) and the arguments 
of the error functions in the integrand. It is convenient 
to analyze the magnetic field response by considering the 
ratio (5I S (B; k))/(SI s (k = 0)), where the denominator is 
given by Eq. 
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where z = B/bi is the scaled magnetic field. For nota- 
tional convenience we introduced the K-depcndcnt terms 
Q and 7, which are defined as 



b 2 TT D = 1 J TO. = baT 



c = 
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2(1 -K 2 ) 



(1 - 

2k 2 (1 + k 2 )' 



7 = 



1-k 2 



1 



1 



(1 - K 2 ) 3 
K 2 (l+K 2 ) 



(74) 



(75) 



It is seen that the arguments of the error-functions as 
well as the power in the exponent diverge in the limit 
k — > 1, i.e. when the <?-factor of one pair-partner is zero. 
This divergence signifies that magnetic field response is 
weak for small (1 — n). The underlying reason for this is 
that the portion of soft pairs goes to zero if the levels of 
one of the partners are not split by a magnetic field. In 
Fig. [7] we plot the magnetic field response for different 
values of k. There are two noteworthy features of this 
response. Firstly, the sign of response is opposite to that 
for inequivalent distributions of electrons and holes, see 
Fig. |51 Secondly, the shape of 5I S (B) is not Lorentzian 
anymore. In fact, this shape is close to Gaussian, as il- 
lustrated in the inset. Another peculiar feature of SI S (B) 
which can be seen from Fig. [7] is that, for k close to 1, 
the response SI(B) develops a bump. 



E. Averaging in the fast-recombination regime 



Turning to Eq. (|46j) for t in the fast-recombination 
regime we notice that the second term in the square 
brackets has exactly the same form as the contribution of 
the soft pairs to 1 in the slow-recombination regime, see 
Eq. (|4"5")l . The underlying reason is that, similarly to soft 
pairs, this second term also comes from the slow eigen- 
mode. The origin of this slow cigenmode, i.e. orthogonal- 
ization of S'-mode to all the other states, was discussed in 
detail in Sect. lie. Since the configurational averaging for 
soft pairs was already carried out, we conclude that the 
magnetic field response in the fast-recombination regime 
is simply described by Eq. (IQGD . 

At this point we note that configurational averaging 
over slow pairs was based on the applicability of the con- 
dition 6qTt d ^> 1. Therefore, it is important that this 
condition is compatible with fast-recombination, b$T <C 
1 , by virtue of a small parameter t/t d . 

In addition to the soft-pair contribution, Eq. (|4"6]) 



nator. This term becomes large when fii and Jl 2 are 
collinear. However, the statistical weight of these con- 
figurations is smaller than the statistical weight of the 
soft-pair contribution. Indeed, in order for the term with 
I Oi x O2 1 2 in denominator to become large, the angle 
between the vectors fii and 1~2 2 should be restricted to 
l/^o \/ tt d "C 1. In course of configurational averag- 



emerges which is small as 



ing, the integral, J eft? sin # . . 
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FIG. 8: (Color online). Different domains on the plane 
(boT D , bar) illustrate the regions where different OMAR 
mechanisms dominate. The is no OMAR in the white do- 
mains. The pink domain corresponds to slow recombination, 
and OMAR is given by Eq. (|52|) . In both the upper and the 
lower parts of the gray domain the OMAR is dominated by 
soft pairs and is described by Eq. (|66[) . The green line divides 
the gray domain into subregions where the recombination is 
slow (upper part) and fast (lower part). The boundaries of 
the domains are: b r = , 1 , and boT = (tore) 1 ' 3 . 



We now turn to the limit of very weak hyperfine fields 
for which the parameter 6qTt d is small. One may expect 
that magnetic field response is suppressed in this domain. 
What we demonstrate below is that this suppression is 
anomalously strong. Namely, the first term of the expan- 
sion of Eq. (|46|) with respect to b\rT D docs not contain 
the external field at all. This first term has the form 



TTt, 

t-2r D = 

16 



|Oi| 2 - |n 2 | 2 ) 2 , 4|n x xfi 2 | 



|fti + tt 2 | 2 



|fii + tt 2 | 2 



(76) 

To realize that B drops out of the expression in the square 
brackets it is convenient to first replace 1 1 x 1 2 by 
I S~ii 1 2 1 X~22 1 2 — (fii ■ J7 2 ) 2 and then use the identity 

|n 1 +S7 2 | 2 |r2 1 -S7 2 | 2 = (|n 1 | 2 + |f2 2 | 2 ) 2 -4(O r ^ 2 ) 2 . (77) 

This leads to a drastic simplification of Eq. (f76")) . which 
assumes the form 



t-2T D = —£\n 1 -n 2 \ 2 . 



(78) 



Since |Oi — Oa| = | be — &fc|> the magnetic field drops out 
of t in the first order in TT D b\. 



VI. CONCLUDING REMARKS 

(i) Our findings can be summarized in the form of do- 
mains on the plane (boT Dl boT), as shown in Fig. [8] 
The fact that for small 6 2 tt d the OMAR response 
is absent is reflected in Fig. [5] by leaving the domain 
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lying below the hyperbola uncolorcd. Large hyper- 
fine fields, &ot > 1, correspond to slow recombina- 
tion. As we have demonstrated above, the OMAR 
for boT > 1 can be dominated either by "typical" 
pairs or by "soft" pairs. The corresponding regions, 
I and II, are colored in Fig. |8]by pink and gray, re- 
spectively. The domains are separated by the curve 
b r = Oot d ) 1/3 . Eq. (JSll) describes OMAR in the 
domain /, while in the domain 77 Eq. (|66[) applies. 
Note that in the domain II only the part above the 
green line corresponds to slow recombination. The 
part below the green line corresponds to fast re- 
combination, but Eq. (|66|) applies in both domains. 
The diagram describes the regimes of OMAR in low 
applied fields, B ~ bo- As B increases above bo, the 
gray domain shrinks. 

(ii) The OMAR response from the soft pairs relies ex- 
clusively on the asymmetry between electron and 
hole. The evidence in favor of such an asymme- 
try was inferred in Ref. from the analysis of 
magnetic-resonance data in organic devices. In Ref. 
l27l . the ratio b 2 /bi was estimated to be close to 3, 
which leads to the value of the asymmetry param- 
eter rj « 0.9. Note, that bipolaron mechanism is 
insensitive to the asymmetry between electron and 
hole. 

(iii) "Parallel-antiparallel" mechanism of Ref. UU yields 
the OMAR response on the level of rate equations 
with the transition rates calculated from the golden 
rule. The applicability of this treatment requires 
that the separation of Zceman levels is large com- 
pared to their widths. On the other hand, the 
OMAR response based on soft pairs, studied in the 
present paper, comes entirely from pairs for which 
the Zeeman levels are almost aligned. This requires 
one to go beyond the golden rule. Previously, a sim- 
ilar situation was encountered 2 ^ by M. Schultz and 
F. von Oppcn in the study of transport through a 
nanostructure with almost degenerate levels. The 
role of spin-selective recombination was played by 
coupling to the leads which was strongly different 
for symmetric and antisymmetric combinations of 
the wave functions. M. Schultz and F. von Oppcn 
pointed out that when two levels are closer in en- 
ergy than the width of each of them, then the con- 
ventional rate-equation-based description is insuffi- 
cient. 

On the physical level, the near-degeneracy implies 
that some spin configuration is preserved during 
many precession periods, i.e. the dynamics is im- 
portant. To account for dynamics, it is intu- 
itively appealing to take the result of Schulten and 
Wolynes, Eqs. ([2])-® , and multiply it by a fac- 
tor describing exponential decay of population of 
states due to recombination. Such an approach was 
adopted in Ref. H^. What this approach misses is 



the feedback of recombination on the pair dynam- 
ics. It is the central message of the present paper 
that this effect is strong in certain regimes, since 
feedback creates long-living modes. 

(iv) The "parallel-antiparallel" mechanism of Ref. UH is 
based on the picture of incoherent hopping of one 
of the charge carriers on the site already occupied 
by the other carrier. We considered the transport 
model applicable for bipolar system where the pas- 
sage of current is due to recombination of electrons 
and holes. However, the principal ingredients of 
both models are the same: (a) in both transport 
models the spins of the carriers precess in their effec- 
tive magnetic fields, the precession being governed 
by the same Hamiltonian Eq. ([I]); (b) the passage 
of current is the sequence of cycles, only one step 
of each cycle is sensitive to the spin precession; (c) 
whether it is a hop or recombination, it occurs only 
from the S'-spin configuration; (d) if either the hop 
or recombination act takes too long, the carriers by- 
pass each other. 

(v) Both the "parallel-antiparallel" pairs and soft pairs 
create the OMAR response by blocking the current. 
The origin of this blocking is completely different 
for the two mechanisms. In the former, the cur- 
rent is blocked due to collincarity of full fields for 
the pair-partners, while for the latter the block- 
ing is due to coincidence of their absolute values. 
In general, both contributions are present in the 
fast-recombination regime. The contribution of soft 
pairs in this regime dominates by virtue of their 
statistical weight. 

(vi) Another distinctive feature of the soft-pairs mech- 
anism follows from Eq. ([55)1 . It contains a combi- 
nation (|fii| 2 — |f22| 2 ) 2 in the denominator. As the 
precession frequencies change with external field, B, 
the pair undergoes evolution from typical to soft 
(when I = |fi2|) and back to typical. Impor- 
tantly, this evolution takes place within a narrow 
interval of B, so that at a given B only certain 
sparse pairs contribute to the current. As demon- 
strated in Ref. [H, this redistribution of soft pairs 
gives rise to mesoscopic features in 1(B) in small 
samples. 

(vii) We have demonstrated above that regardless of 
whether the OMAR is due to blocking caused by 
"parallel-antiparallel" configurations, as in Ref. [jjl, 
or due to soft pairs, the shape of the response is 
always close to B 2 /(B 2 + B 2 ). This result was 
obtained under the assumption that r and t d are 
fixed. If the values of t and t d are broadly dis- 
tributed, then the adequate description of trans- 
port should be based on the percolative approach^ 2 .. 
However, within our minimal model, the current is 
the sum of partial currents through the chains, see 
Fig. [3] Then, with wide spread in cycle durations, 
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t, the current will be limited by pairs with longest 
t present in each chain. 
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with dissipation term redefined as L(p) = 
-\T{p, \S) (S\}. To derive Eq. dHJ), we search 



(A9) 



for solution of Eq. (|A8j) in the form 



and find that ?/>(i) must satisfy the following 11011- 
hermitian Schrodingcr equation 



i^\m) = H'\m), 



(A10) 



Appendix A: Time Evolution and the Schrodinger 
Equation 

In this Appendix we sketch a formal derivation of Eqs. 
(TT9"]) and ([57)1 starting from the Liouville equation for the 
density operator, <r, 



da_ 



-i[H,a] +L{a), 



(Al) 



where the term L(a) describes relaxation, which in our 
case is recombination from S to the ground state, G. 
The ground state with energy —£ is included into the 
bare Hamiltonian 

H = H + H G (A2) 
= (s l -B 1 + S 2 -B 2 ) -£\G)(G\. (A3) 



Then the operator L(a) cast into conventional Lindblad 
form2£ reads 

L{d) = \T (2 \G) (S\ a \S) (G\ a \S) (S\ \S) (S\ a) , 

(A4) 

where = r 1 is the inverse recombination time. 

Denote with i, k different spin configurations of the 
pair prior to recombination. The form Eq. (|A4[) of the 
dissipation ensures independence of the elements of the 
density matrix with subindices z, k from the elements 
containing subindex g. This decoupling follows from the 
full system of the equations of motion 



<9ov 



dt 
da Gk 



dt 
da, 



rcr ss , 



(A5) 



- -jFoaJsk, (A6) 



= -i[H,a] ik --T{a,\S) (S\} ik 



(A7) 



Eq. (|A7j) couples only the elements of 4 x 4 matrix, which 
we denote with p, so that Eq. (|A7I) represents equation 
of motion for p. These equations can be rewritten in the 
form similar to Eq. (|A1|) 



dp 
~dt 



-i[H,p\+L(p), 



(A8) 



where H' is defined as H' = H - if \S) (S\. The fact 
that decoupling Eq. (|A9|) is valid follows from a straight- 
forward calculation 



4 w m = w) m + iv>> (i^- t 



= [H-i-\S) (VI 



|V) ty\[H + i-\S) (S\ 



(All) 
(A12) 

(A13) 



= [H,W)m-i-{\S)(S\,W)W}. (A14) 

Now Eq. (|19[) immediately emerges as an equation for 
eigenvalues of the operator H' . 



Appendix B: Derivation of Eq. (|37|) 

To derive Eq. (|37[) for recombination time from ran- 
dom initial state, we first find the expression for recombi- 
nation time, t^ , from a given initial state, ipo , in terms of 
the solution of Eq. (|A8|) for p(t) complemented with con- 
dition p(0) = \ipo) (Vol- The expression for t^ in terms 
of the full density matrix a(t) reads 



t 



(Bl) 



The meaning of the expression in the brackets is the prob- 
ability that recombination took place between t and t+dt. 
The expression for t^ in terms of p(t) follows from the 
relation 



°~gg + Trp = 1. 
Performing integration by parts, we obtain 

/>oo 

t^ = / dtTrp(t), 



(B2) 



(B3) 



To find the recombination time (t n ) from the random 
initial state the time t^ should be averaged over initial 
states. One way to perform this averaging is to fix a 
certain orthonormal basis, expand ipo as 



Vo = ^c fe $ fe , 
k 



(B4) 
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and express tj, Q as bilinear form in Cfe. This yields 



1 4>Q 



dt Tr 



diTr 



(B5) 



CT(t) 



E c fe' $ fe' 
. fe' 



(B6) 
(B7) 



where {/ (t) is the non-unitary evolution operator. Now 
the averaging over initial conditions reduces to averag- 
ing over Cfc according to the rule (cfcC$) = \Sk,k'- This 
averaging is straightforward leading to 



(tn) 



dt Tr 



(B8) 



The remaining task is to express the sum, Eq. (|B8[) . in 
terms of eigenvalues and eigenvectors of a non-hermitian 
Schrodinger equation, Eq. (|A10|) . To accomplish this 
task we will use the expansion of the solutions ip k of Eq. 
(|A10[) . which we, for brevity, denote with |A fc ), in terms 
of the orthonormal basis $ fc , which we denote with \k). 

In terms of these new notations Eq. (|A10|) and the 
time evolution operator can be written as 

H' \Xj) = Xj , U(t) \Xj) = e- 1 ^ \Xj) . (B9) 

It is also convenient to introduce a matrix, d, which 
relates the elements of the basis to the solutions of Eq. 
(|AlQ|) . Namely, 

\k) =J2dkl\Xi) . (BIO) 
i 

Substituting Eq. ([BT0|) into Eq. (HU), we find 

U(t) \k) = d k .iU(t) |A») = £ d kl e^ 1 |A,} . (Bll) 
( i 

Next we introduce, g 7 which is the matrix of scalar prod- 
ucts 

g ij = (X i \X j ). (B12) 
Using the definitions Eq. (|BT0|) and Eq. (fBT2l) we express 
(t R ), defined by Eq. (|B8|) . in terms of the matrices d and 



(**) 



5 E J dtTv[U(t)\k) (k\uHt) 

k g 

iE/°° rfiTr ^E^i A *>) 



(B13) 



, (B14) 



OO 

JE / dte-^- x '^d kl dl m Tr[\X l ){X m \] 1 



(B15) 
(B16) 



In the last identity we have isolated the combination of 
the elements of the matrix d. The reason is that this 
combination can be cast in the form 



^2d k id* km - g ml 



(B17) 



To prove the latter identity, we start from the matrix 
relation 



J Xi I i) = E dij ( X i I A i> = E di o9lj > 



(B18) 



and invert it to obtain 



(B19) 



Next we complex conjugate both sides of Eq. (|B19j) 
which yields 



4 = E-^ 1 * < A <i 0* - E-^ 1 * (*i A <> • ( B2 °) 



Now, the identity Eq. (|B17[) emerges as a result of 
straight-forward calculation 



k k \ j J \ n / 

(B21) 

= E%V7™ (Xj\ ( E l fc ) < fc l ) l A ™> (B22) 
= E% 1 9„„*3in (B23) 



-1* 



9ml 



(B24) 



Finally, substituting Eq. (|B17|) into Eq. (|B16|) , we arrive 
at Eq. (|3T|) of the main text . 
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